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Abstract. 

The electron Bernstein wave (EBW) is typically the only wave in the electron 
cyclotron (EC) range that can be applied in spherical tokamaks for heating and current 
drive (H&CD). Spherical tokamaks (STs) operate generally in high-/? regimes, in which 
the usual EC O- and X- modes are cut-off. In this case, EBWs seem to be the only 
option that can provide features similar to the EC waves — controllable localized H&CD 
that can be utilized for core plasma heating as well as for accurate plasma stabilization. 

The EBW is a quasi-electrostatic wave that can be excited by mode conversion from 
a suitably launched O- or X-mode; its propagation further inside the plasma is strongly 
influenced by the plasma parameters. These rather awkward properties make its 
application somewhat more difficult. In this paper we perform an extensive numerical 
study of EBW H&CD performance in four typical ST plasmas (NSTX L- and H-mode, 
MAST Upgrade, NHTX). Coupled ray-tracing (AMR) and Fokker-Planck (LUKE) 
codes are employed to simulate EBWs of varying frequencies and launch conditions, 
which are the fundamental EBW parameters that can be chosen and controlled. Our 
results indicate that an efficient and universal EBW H&CD system is indeed viable. In 
particular, power can be deposited and current reasonably efficiently driven across the 
whole plasma radius. Such a system could be controlled by a suitably chosen launching 
antenna vertical position and would also be sufficiently robust. 
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1. Introduction 

Present research in electron cyclotron (EC) wave heating and current drive (H&CD) for 
magnetic confinement thermonuclear fusion [1] is focused on conventional aspect ratio 
tokamaks, and particularly ITER. However, in the "alternative" spherical tokamak (ST) 
with aspect ratio A = Ro/a close to unity {Rq and a being the major and minor radii, 
respectively) and weaker external toroidal magnetic field, the usual EC transverse O- 
and X-modes are mostly cut-off and cannot be used for H&CD. The role of ST research is 
nevertheless very important. For their relatively high neutron fiux density and economy, 
STs are being considered as a candidate for a component test facility (ST-CTF) [2], |3] 
and, for the same reasons, appear in fusion-fission hybrid concepts [4j. 

ST's low magnetic field has a major impact on the propagation of electron cyclotron 
waves in the plasma. This frequency range is of crucial importance for auxiliary 
H&CD systems in present and future tokamaks. Typically, in STs, the electron plasma 

1 /2 

frequency oope = 27r/pe = (neC^/meeo) is much greater than the electron cyclotron 
frequency Uce = 27r/ce = eB/nic- Here, is the electron density, B is the total magnetic 
field, e is the electron charge, rrie is the electron mass, and Eq is the vacuum permittivity. 
A similar situation often arises in stellarators, which do not have any principal MHD 
stabihty density limits. In this so called overdense regime, particularly when cUpe > nuce, 
n > 3 in most of the plasma cross section, the O- and X-modes of EC waves with 1*^*, 2'^'^ 
and 3'''^ harmonic frequencies (wce < OJ < 4:0Jce) are cut-off and cannot propagate inside 
the overdense plasma. Higher harmonic EC waves are not of interest because of their 
very weak absorption (low optical depth). However, the electron Bernstein wave (EBW) 
[5] — a quasi-electrostatic kinetic EC mode — can propagate and be strongly absorbed in 
an overdense plasma. 

EC waves are extremely useful because they can be launched far from the 
plasma (they do not need large plasma-facing antenna structures like ion-cyclotron or 
lower-hybrid waves) and feature highly localized and controllable H&CD capabilities. 
The application of the overdense mode — the EBW — is, however, complicated by its 
electrostatic nature. First, EBWs must be excited by appropriately launched O- or 
X-mode via so called OXB or XB mode conversion scheme. This mode conversion 
takes place in the upper hybrid resonance region, where the wave frequency satisfies 
u = = a/i^pc + ^ce- This typically occurs near the plasma edge. The mode 

conversion efficiency depends on the wave and plasma parameters and is thus a potential 
source of power loss. The excited EBW can subsequently propagate inside the overdense 
plasma; however, because of its dispersion characteristics, the propagation strongly 
depends on plasma parameters and the wave vector can change considerably in various 
ways (unlike O- and X-mode propagation, during which the parallel wave number is 
mostly conserved). 

H&CD by EBWs have already been demonstrated experimentally in magnetic 
confinement fusion devices, particularly in COMPASS-D [6] and Wendelstein 7-AS [HIH]. 
Numerical studies of advanced (steady-state) spherical tokamak operations consider 
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EBWs as one of the vital current drive system that can stabihze MHD instabihties (H]. 

In this paper, we pursue an overall study of EBW H&CD on spherical tokamaks 
by means of numerical ray-tracing and Fokker-Planck simulations. Two coupled 
codes— AMR (Antenna, Mode- conversion, Ray-tracing) [TOl [11] and LUKE [121 113]— 
are employed. These codes have proven to be very suitable for EBWs [HI [TT]. A large 
number of cases with different injection parameters is simulated in four different ST 
conditions: two experimental discharges of the NSTX tokamak [15], an ST-CTF-like 
MAST-Upgrade H-mode scenario [HI [IT] and an NHTX scenario [15]. The resulting 
extensive collection of results is analyzed with emphasis on H&CD performance — 
viability, effectiveness, flexibility, controllability and robustness. 

The paper is organized into five sections. After the introduction, the description 
of EBW physics involved in this study is given. Section [3] describes the context of 
the numerical simulation — the target plasmas and the EBW parameters. Numerical 
results — particularly the current drive localization and efficiency, the role of N\\, Z^fi, 
the quasilinear effects and the robustness — are presented in section HI Finally, section 
[5] discusses the conclusions of our work. 

2. Electron Bernstein wave description 

The electron Bernstein wave (EBW) has been known for decades [5j. It appears as 
another EC branch, besides the cold plasma O- and X-modes, when solving the kinetic 
dispersion relation of a plasma in an external magnetic field. It was later applied 
to magnetic confinement fusion for H&CD while EBW emission was introduced as a 
diagnostic tool; an overview of EBW experiments is given in [18]. EBW physics can 
basically be separated into three areas: the mode conversion, the propagation and the 
wave-plasma power transfer. 

2.1. Mode-conversion 

The EBW mode conversion is a full wave process in which transverse electron cyclotron 
modes and the EBW are involved. The mode conversion always requires the slow X- 
mode to be excited, which can then fully convert into the EBW at the upper hybrid 
resonance (UHR). This study is confined to EBW excitation from the low- field side 
(LFS) of a tokamak, where two possibilities exist: the XB [19j and the OXB [201 [21] 
conversion. The XB conversion is characterized by a direct coupling between the fast 
and the slow X-mode while in the OXB scheme the 0-mode is converted to the slow 
X-mode (in fact, it is converted to the fast branch of the X-mode, which propagates 
towards higher density and then smoothly converts to the slow branch, which propagates 
backwards to the UHR). The XB scheme is typically efficient for lower frequencies and 
requires the density scale length to be specifically adjusted (for details see [IH]). The 
OXB scheme, on the other hand, is more universal in terms of frequencies and density 
scale lengths as efficient conversion only requires the 0-mode to be incident at the 
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optimum angle. In ID theory, this angle is given by [201 123] 



(1) 



where N = ck/u; is the normalized wave vector, A*"]! = N ■ B/i?, and A'^poi = 
N ■ (B X V^c) / ||B X Vne||. If the incident angle is not optimum, the O- to X-mode 
power conversion efficiency of a plane wave decays approximately exponentially as 



where L„ = Ue/ (dne/dr) is the density scale length, and (j2]) is evaluated at the 0-mode 
cut-off. This formula agrees well with numerical results for koLn > 1 when N\\ is almost 
optimum [22] so that Cqx > 0.9. When N\\ is further from the optimum, larger koL^ is 
required for good agreement [22]. A ID formalism has also been developed [21] which 
solves the mode-conversion problem for almost arbitrary inhomogeneity scale length. 
In general conditions, the conversion efficiency must be calculated by full-wave codes, 
either in the cold plasma [25] or the hot plasma [26] approximation. These codes also 
take into account the incident wave polarization and the parasitic slow to fast X-mode 
tunneling, which can decrease the OXB conversion efficiency, especially when the density 
gradient is very steep and consequently k^Ln is small. Numerical and analytical results 
are in good agreement particularly around the optimum incidence. As a consequence 
of the exponential efficiency dependence, there always exists an angular window where 
the mode conversion is sufficiently effective. 

Recently, 2D theory and simulations of the OXB conversion have been developed 
[23 1211 1211 [30]. 2D effects are shown to be important for off-equatorial launch, where 
the 0-mode cutoff and X-mode cutoff surfaces are no longer parallel [29]. In [30], it is 
shown that the beam curvature should be matched to the plasma surface curvature in 
order to optimize the conversion efficiency. The effect of the beam size, which obviously 
determines the beam spectrum, is also studied, showing that larger beams, i.e. narrower 
fc-spectrum, tend to be more efficiently converted [30]. Non-linear effects can also play 
a role, for example a parametric decay [31] or higher harmonic wave generation [32] . 
Another factor typically not considered in the conversion process is density fluctuations. 
An estimate, based on a probability distribution function and the ID formula ([2]), was 
made and experimentally demonstrated in [^, showing a considerable decrease of the 
conversion efficiency for large density scale lengths. A detailed treatment of this problem 
should employ a 2D approach. 

In this paper, even though the mode conversion is not treated in detail, it is 
nonetheless not neglected. We consider the OXB scheme for its universality. However, 
the results are directly applicable to any mode conversion scheme as long as EBWs 
are excited at the same place with similar N\\. The OXB scheme was successfully 
demonstrated in various past and present experiments [l8j. Our EBW H&CD simulation 
starts from an antenna, which emits a Gaussian beam [33] of a given frequency and waist 
radius wq, which, in our case, is calculated from the Rayleigh range 



[211 123] 




(2) 



zr = Tcwl/Xo = kowl/2 



(3) 
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where Aq and ko are the vacuum wavelength and the wavenumber, respectively. At the 
distance from the waist, the Gaussian beam doubles its spot size (the beam radius 
becomes \/2wo). At a fixed zr, the beam divergence is similar for all frequencies. Large 
divergence, i.e., wide beam fc-spectrum, would cause poor 0-X conversion efficiency. 

We now calculate the conversion efficiency of a Gaussian beam. The electric field 
of a Gaussian beam in the Fourier space is 

E^exp(^-^{kl + kl)y (4) 

where k^^y are wave vectors perpendicular to the direction of the beam propagation. 
The corresponding energy density is oc E . We can evaluate the total power of the 
0-X converted Gaussian beam using the analytic formula and Parseval's theorem: 

Pox ^ /■!l/.!ll-g|'CoxdMfci; 

Po r r lEfdk^dky ■ ^ ^ 

We assume here that the error introduced by changing the integration limits from 
k'^ + ky ^ /cq to ±oo is negligible, which is valid for beams with not too large divergence. 
We now assume that the z axis is along the beam propagation direction and that the 
launch angle is optimum by setting, without any loss of generality, the magnetic field at 
the 0-mode cut-off to be B/i? = ^0, ±y^l — ^l^opt' ^ll.opt) • ^-^-^ magnetic field is 
assumed to be homogeneous across the beam spot at the 0-mode cutoff. The integral 
(El) then becomes 



ox 



J —OO J —oc II X \ 

J — oo J — oo II y 

This can be evaluated analytically after Taylor-expanding the exponent in kx^y around 
to second order, yielding 

Pox .ft:'^"^ 



where k = 7rLn\/uj^^j2uj. This is an important result which, in fact, imposes an upper 
limit to the conversion efficiency of a Gaussian beam. This limit depends on L„/ zr for a 
fixed ooce/oo. (Note that 0.4 < a/wcc/2u; < 0.7 for the first two harmonics). It also tells us 
how narrow (i.e. how divergent) a beam can be used while keeping the OXB conversion 
efficient. The beam conversion efficiency is shown graphically for the discussed scenarios 
in section 3.2. In a similar fashion, we can also evaluate the conversion efficiency of a 
Gaussian beam for non-optimum central wave vector. 



2.2. Propagation 

The electron Bernstein wave is, apart from the "cold" O- and X-modes, a solution 
to the kinetic (hot) dispersion relation of a plasma in an external magnetic field [5]. 
Numerous analytical and numerical studies of EBW propagation have been performed 
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and EBW propagation is hence quite well explored. The characteristic properties of 
EBW propagation are: 

• The polarization is quasi-electrostatic; hence, in most situations, the electrostatic 
dispersion relation describes EBWs satisfactorily [M f |35 | [36]. 

• EBWs can propagate in plasmas if lOcc < ^ ^ i^uiJil- As w ~ nuce, n = 1,2,..., it 
implies nUce < c^uh- Note again that EC O- and X-mode are mostly cut-off under 
these conditions. There is no upper density cut-off for EBWs. 

• The phase velocity is almost perpendicular to the external magnetic field {k± ^ 
the group velocity is generally different from the phase velocity in both magnitude 
and direction. 

• The wavelength is of the order of the electron gyroradius pe = v^e/uce {v^e = 
^jTf.jm^ is the electron thermal velocity and is the electron temperature in 
energy units), i.e., 

fcPe = fc±Pe ~ 1. (8) 

• EBW characteristics vary significantly depends on whether it approaches a 
resonance at the lower harmonic (cj > nWce) or at the higher harmonic (w < nWce)- 
In particular, the perpendicular wave vector is much smaller near the resonance 
if cu < ncjce, where k^p^ <C 1 and electromagnetic effects are no longer negligible. 
Yet, as the wave approaches the resonance, the power is usually absorbed before 
the electrostatic approximation becomes invalid [Mj. 

• The parallel refractive index iVy evolves during the propagation and can be greater 
than one (unlike O- and X- modes). Depending on the magnetic field topology 
and the vertical launch position, the wave parallel index can either stay close to 
its initial value, or oscillate around zero, or increase/decrease steadily. During 
the actual EBW propagation, the wave parallel index can change the regime. In 
general, waves close to the midplane tend to have a flat or oscillating N\\ [37] , while 
for off-midplane rays N\\ increases or decreases steadily at a rate proportional to 
the distance from the midplane [19]. This property forms the basis for controlling 
EBWs. 

EBW propagation in a tokamak plasma is far from trivial and necessitates a numerical 
simulation. The ray-tracing technique is well suited for EBW propagation since the 
WKB validity conditions are well fulfilled due to the short wavelength. We employ the 
AMR code [101 E] to simulate the EBW propagation. This code uses a conventional ray- 
tracing method |38l [39] with an electrostatic kinetic non-relativistic dispersion relation 

m 

I w < wuH holds for u) < 2llIco- Higher frequency EBWs can propagate in a rarified plasma up to a 
certain lower electron density limit. 
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= nrf, i , ^ = , (10) 



where 



Z is the plasma dispersion function |?T] and J„ is the modified Bessel function of the 
first kind. The ray trajectory is then a solution to the Hamiltonian-type equations 
dr _ _d^{V) d^{V) 

dt dk du ^ /-riN 

dk _ d^jV) d^(V) ^ ' 

dt dr duj 

Equations (ITTi) conserve I' = at zeroth order in the parameter S3 (P) /3fJ {T>), i.e., in a 
weak damping approximation. 

In this study, we neglect the collisional damping, which can, however, be of 
critical importance, as was previously shown by our modelling and by experiments at 
NSTX [IJlHSj. Collisional damping is extremely sensitive to edge plasma conditions, 
particularly the temperature (or, more precisely, the collision frequency) in the mode 
conversion layer. As these conditions cannot be sufficiently accurately predicted, we 
completely ignore the effect of collisions. 

By using the non- relativistic electrostatic dispersion ([9]), we neglect relativistic and 
transverse electric field effects in the wave propagation. These effects were studied in 
[miSlj, and also in [45j, where a fuUy-relativistic electromagnetic ray-tracing code was 
used. Some differences are seen in the ray propagation and in the polarization. These 
differences are, however, rather small and would not change the overall picture of our 
survey. The fuUy-relativistic approach is also computationally very intensive and would 
require an extensive amount of computation time to simulate all the cases presented 
here. 



2.3. Wave absorption 

In this section, we describe the theory of the EBW absorption. EBWs can be absorbed 
by resonant electrons, which satisfy the resonance condition 

u — nUcc/j — k\\V\\ = 0, (12) 

— 1 /2 

where 7 = (1 — f^/c^) is the usual relativistic factor. While the absorption 
mechanism — the EC harmonic damping — is identical to that for O- and X-modes, 
the polarization of EBWs is different (quasi-electrostatic), and this leads to a strong 
interaction even for low temperatures at any EC harmonic. In the ray-tracing approach, 
the wave absorption is assumed to be weak and can thus be handled independently. This 
requires the anti-Hermitian part of the dielectric tensor to be much smaller than the 
Hermitian part, or, alternatively, 53 {V) <^ 3? (D). 

The zeroth order solution of the dispersion relation leads, as shown in the previous 
section, to the ray-tracing equations (fTTj) . The first order solution then leads to the 
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radiative transfer equation [3H] 




(13) 



where 



(14) 



\d^{V) /du 



is the absorption coefficient, i!} the emissivity and P the ray power. The non-relativistic 
dispersion function is, however, inappropriate for the damping calculations. For this 
reason, the ray-tracing code rather employs the weakly-relativistic absorption coefficient 
by Decker and Ram [34], which is very fast and sufficiently accurate. Assuming high 
incident ray power Pq, the emissivity can be neglected and the solution to Equation (fT3ll 
is 



Since equation ([T5|) provides a linear damping solution, any modifications to the electron 
distribution function by the waves are not taken into account. However, quasilinear 
effects due to the modified distribution function can play an important role in EC 
H&CD. For this reason, we employ LUKE — a fully relativistic, bounce- averaged, 3-D 
Fokker-Planck solver [12] — which calculates the evolution of the electron distribution 
function for axisymmetric plasmas in the low-coUisionality regime. LUKE particularly 
accounts for collisions and quasilinear diffusion due to RF waves. The code uses a fuUy- 
relativistic, parallel momentum conserving collision operator and a fully implicit 3-D 
time evolution scheme for a fast convergence to the time-asymptotic solution. It has 
been verified that the damping proffie calculated by LUKE in the low power limit agrees 
with linear theory [T3] . 

2.4- Simulation model 

Our model antenna emits a Gaussian beam, parameterized by its frequency, beam waist 
vertical and radial position Za, -Ra and waist radius Wq. The beam radius is, in our 
case, calculated from the Rayleigh range, therefore fixing the beam divergence and 
consequently the beam 0-X conversion efficiency (putting aside the variable density 
scale lengths and the Uce/oj dependence) for all frequencies. 

The antenna angles used in the following simulations are optimized for the OXB 
mode conversion, i.e., determined by the condition ([1]) for each beam waist position. 
The average over a single harmonic frequency range (nwce < uj < {n + 1) ujce) is chosen 
as the optimum angle for these frequencies rather than the optimum evaluated for a 
particular frequency. The differences in the resulting ray-tracing initial conditions are 
negligible and central ray mode-conversion efficiencies do not drop below 90 % in most 
cases. The launch conditions are also checked using AMR's ID full-wave calculation 



The beam is discretized by a bundle of 16 individual rays. We denote the positions 
of the intersections of the rays with the beam waist by r^, z = 1 . . . 16. Using a straight 



(15) 



[25j. 
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line propagation, the intersection of the central ray with the 0-mode cut-off is found. At 
this point, the Gaussian beam size is calculated and the beam is again discretized by the 
same, proportionally scaled 16 ray pattern. We denote these ray spot positions rj. The 
intersections of the straight lines, connecting and r^^, with the 0-mode cut-off surface 
are used as the starting points for the ray-tracing: i.e., we assume straight propagation 
of the 0-mode from the beam waist to the 0-mode cutoff surface, but still taking the 
beam divergence into account. The initial wave vectors for the ray-tracing are found by 
solving the electrostatic dispersion relation with respect to N±, using 

iV||o = No ■ B/B, (16) 

where Nq = (r^^^ — r^) / Hr^^ — r^|| is the ray vacuum wave vector and B is evaluated at 
the 0-mode cut-off. 

After finding the initial wave vector, the electrostatic ray-tracing is started. The 
principal results are the ray trajectories and wave vector evolutions. When ray-tracing 
is finished, the outputs are passed to LUKE, along with the magnetic equilibrium and 
plasma profiles. The AMR-LUKE interface has been particularly verified to keep all 
quantities consistent. Besides, the interface is user friendly and LUKE can be launched 
by AMR, and vice versa, by a single option in the configuration file. 

Finally, LUKE determines a distribution function /qi that is consistent with the 
quasilinear wave absorption. Besides the flux-averaged absorbed power density profile 
Pd (p), the fiux-averaged EBW-driven current density profile 



J\\ (P) = "^(/ ^ll/qidp^') 



(17) 



is also calculated. Here, (■) denotes flux surface averaging and p is a flux surface 
coordinate based on the poloidal magnetic flux ip: 



P = Ppol 



V'LCFS - V^axis 

where ipa.xis and '?/'lcfs is the poloidal magnetic flux at the magnetic axis and at the last 
closed flux surface (LCFS), respectively. 



3. Simulated scenarios of EBW H&CD 

3.1. Fundamental target plasma parameters 

As already mentioned in the introduction, EBW H&CD is simulated here in four 
different target spherical tokamak scenarios, whose fundamental parameters are listed in 
Table [H As can be seen, the chosen scenarios differ in various fundamental parameters. 
Two of them are typical NSTX L- and H-mode experimental discharges, the other two 
are TRANSP (a plasma transport code) [56j model scenarios of the planned MAST 
Upgrade [I7j and of NHTX (a potential plasma facing component test facility). 

In Figure [1] are plotted the midplane radial profiles of the characteristic frequencies 
for the target scenarios. The simulated frequency ranges are marked by shaded areas. 
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Table 1. Fundamental parameters of the target ST scenarios: major radius Rq [m], 
minor radius a [m] , central toroidal magnetic field Bq [T] , central electron density rico 
[10^^ m"'^], central temperature Tco [keV], plasma current Ip [MA]. 



Name 


Rq 


a 


Bo 








Origin 


NSTX L-mode 


1.0 


0.52 


0.5 


2.6 


2.9 


0.6 


shot 123435 


NSTX H-mode 


1.0 


0.5 


0.5 


3.9 


1.4 


1.0 


shot 130607 


MAST Upgrade 


0.93 


0.41 


0.78 


3.5 


2.4 


1.2 


TRANSP 


NHTX 


1.2 


0.37 


2.0 


20.0 


5.7 


3.5 


TRANSP 



whose areas delimit the EBW propagation regions, i.e., from the UHR at the edge to 
the cold EC resonance, which is theoretically accessible by iVy = waves only. Clearly, 
the plasma is overdense (/pe > nfce) and the first three EC harmonics are inaccessible 
for O- and X-modes in NSTX and MAST. In NHTX, the third harmonic is more or 
less accessible and the corresponding frequency is compatible with present-day gyrotron 
technology. The first two EC harmonics have been selected for NSTX and MAST, as 
higher harmonics will likely be overlapping because of the Doppler broadening. The 
same applies to NHTX, where, however, only the first harmonic is simulated since the 
second is only marginally overdense and the OXB conversion region occurs in the core 
plasma rather than at the plasma edge. Moreover, the high second harmonic frequency 
(~100 GHz) combined with the relatively long density scale length in the conversion 
regions makes the OXB mode conversion angular window rather narrow. 

3.2. EBW system parameters 

As already shown in numerous previous works (see, e.g., |171[l9]), the propagation path 
and the iVy evolution of EBW rays strongly depends on the vertical launch position. 
N\\ appears in the resonance condition ( |T2i) and can thus influence the wave absorption 
location. Moreover, is an important factor in the EBW current drive. The vertical 
launch position is therefore a crucial parameter of any EBW launcher, which, besides the 
frequency, can be chosen arbitrarily. The toroidal and poloidal angles must be optimized 
for the conversion efficiency and we can only select negative or positive initial N\\o- 
These properties and restrictions dictate the EBW launcher parameters to be scanned 
in our survey: wave frequency, vertical launch position, and the sign of Np. The waves 
propagate far enough from the plasma top and bottom so that we can assume up-down 
symmetry. It has been verified for several cases that the above-midplane launcher power 
and current density profiles are symmetric with the below-midplane launcher profiles 
with the opposite sign of Np. In particular, 

P{p,N\\o,za) = P {p,-N\\o,-za) , .^g^ 
j {p, Ni\o, za) = -j (p, -A^iio, -za) ■ 
The reason for this is that the fiux surfaces are axially symmetric with respect to the 
midplane. As a result, EBW rays below and above the midplane are driven in opposite 
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0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 
R [m] 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 
R[m] 




MAST-U 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 
R [m] 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 
R[m] 



Figure 1. Radial profiles (in the mid-plane) of the characteristic frequencies for the 
target scenarios, fpc - electron plasma frequency, /uhr - upper hybrid frequency, nfce 
- n'^ EC harmonic. Filled areas represent the simulated frequency ranges. 



vertical directions with N\\ having opposite signs. Subsequently, the ray paths are 
symmetric when starting from ±za with opposite N\^, and thus resulting in symmetric 
power deposition profiles. The opposite sign of N^^ then causes a reversal in the driven 
current direction. 

The frequencies and vertical launch positions used in our survey are given in Table 
Table [2j The antenna beam Rayleigh range for all simulations is set to 0.5 m. In Figures 
[2] and [3] are shown the maximum Gaussian beam 0-X conversion efficiencies calculated 
from Equation ([71) for typical target plasma density scale lengths and the dependence on 
Ln for various Rayleigh ranges. The selected zn of 0.5 m is a compromise between small 
beams with poor conversion efficiency and large beams with high conversion efficiency 
but presumably wide power deposition profiles. 

A launcher that would span the whole parameter range listed in Table [2] is certainly 
unrealistic. Multi-frequency systems are rather challenging — even though such systems 
are actively being studied [1^. For these reasons we focus more on single- frequency 
systems. A concept of such a system is sketched in Figure HI This system is designed 
to have a vertically movable mirror, which can be rotated in two dimensions (toroidally 
and poloidally), thus providing variable vertical launch positions with optimum OXB 
launch angles at the same time. In-vessel components are shielded from the plasma by 



A survey of EBW H&CD potential for spherical tokamaks 12 
Table 2. EBW launcher system parameters used in this study. 



scenario 1^* harmonic frequen- 2"^ harmonic frequen- vertical launch posi- 

cies [GHz] cies [GHz] tions [m] 



NSTX L-mode 


11, 


11.5, 12, 


12.5, 


13, 


21, 


22, 


23 


24, 


25, 


26, 


0, 


0.1, 


0.2, 


0.3, 


0.4, 


0.5 




14, 


15, 16, 17, 


18 




27, 


28, 


29, 


30 


















NSTX H-mode 


12, 


12.5, 13, 


13.5, 


14, 


22, 


23, 


24 


25, 


26, 


27, 


0, 


0.1, 


0.2, 


0.3, 


0.4, 


0.5 




14. 


5, 15, 16, 17, 18, 


19, 


28, 


29, 


30 






















20 
































MAST-U 


19, 


20, 21, 22 


, 23, 


24, 


32, 


33, 


34 


35, 


36, 


37, 


0, 


0.1, 


0.2, 


0.3, 


0.4, 


0.5, 




25, 


26, 28 






38, 


39, 


40, 


41 






0. 


6 










NHTX 


52, 


54, 56, 58 


, 60, 


64, 


none 










0, 


0.2, 


0.4, 


0.6, 


0.8, 


1.0, 




68, 


72, 76 


















1. 


2 
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screens that can either be part of the machine vacuum vessel (the smaller vessel variant) 
or placed inside a large MAST-like (cylinder) shape vessel (the larger vessel variant). 
This system is feasible with present day technologies and provides enough flexibility 
required for an advanced EBW system, as will be shown hereinafter. 

poloidal cut toroidal cut 




Figure 4. A possible concept of an EBW launctier system witli a vertically movable 
antenna. The design is sketched for two (existing) vacuum vessel shapes — NSTX-like 
(spherical) and MAST-like (cylindrical). Two variants are proposed — either inside a 
larger vessel or outside a smaller vessel. 

The principle of using EBW control based on a vertically adjustable launcher is 
demonstrated in Figure [5l It is clearly seen that the vertical launch position (besides 
the frequency) has a major effect on the ray propagation and strongly influences the 
location of the power deposition. We also notice the typical behaviour of EBW rays: 
Midplane rays with frequencies for which the cold EC resonance surfaces (typically 
concave shaped) occur in the inboard half of the plasma typically propagate straight until 
they reach the vicinity of a cold EC resonance, where their |iV|| | grows exponentially and 
the rays are damped |37]. Rays with lower frequencies, whose cold EC resonance layers 
appear in the outboard half with typically convex shape, oscillate around the midplane 
and their N\\ oscillate around zero [37]. Rays launched off- midplane are characterized 
by steadily and monotonically varying N\\. This behaviour is shown graphically in 
Section 14.21 This results in a significant Doppler shift of the EC frequency and hence 
these waves are absorbed quite far from the cold resonance. 

4. EBW H&CD performance — numerical results 

In this section, we present our result on EBW H&CD performance in the previously 
listed scenarios. The mode conversion, which is a well-separable problem, is always 
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Figure 5. Ray trajectories for various frequencies and vertical launch positions for 
a) NSTX L-mode, harmonic, b) NSTX L-mode 2"^^ harmonic, c) MAST-U 1"* 
harmonic. Ray trajectories are plotted with solid lines, dashed and dash-dot lines 
show I''' and 2"*^ cold EC resonance surfaces, respectively. 



assumed to be 100%. There are two fundamental reasons for this assumption. First, 
we would like to focus on EBW propagation, absorption and current drive. These 
aspects, as will be shown hereinafter, are complex enough even when neglecting the 
mode conversion. Second, the theoretical description of the mode conversion process is 
so complex that we can arrive to almost any results depending on which theories and 
processes are considered. Moreover, the knowledge of the edge plasma, where the mode 
conversion occurs, is very limited. 

4.I. Localization and current drive efficiency 

The current drive efficiency can be expressed in various ways. Most straightforward and 
suitable from the experimental and engineering point of view is the absolute efficiency 

V ^ ^, (20) 

where Jrf is the total current driven by the RF waves and Pq is the total injected RF 
power. However, the current drive efficiency unavoidably depends on plasma parameters, 
particularly the coUisionality, and hence a quantity that reflects this intrinsic behaviour 
would be better suited for comparing different plasmas and current drive mechanisms. 
Commonly used for EC waves is the normalized efficiency ( |49], which scales out the 
electron density and temperature coUisional effects and the plasma size. However, ( 
does not reflect the intrinsic effects of particle trapping and effective ion charge (Zes). 
The original deflnition assumes that the power deposition proflle is well localized so the 
plasma parameters (density and temperature) do not change there. This is not always 
valid and we therefore use an absorbed power weighted average: 

e^Ro f (p) dJRF ^ 
C = —-^ / / X . dp. 21 
^0 ^0 J P dp 
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Here, cI/rf is the RF driven current in a plasma surface enclosed by p and p + dp flux 
surfaces. In the numerical simulations, LUKE selects a p-grid (based on the power 
deposition proflle) so the volumes become flnite: 

T, ^ V f^^^^, ^^±1^V (22) 



where V (p, p') denotes the plasma volume enclosed by the flux surfaces p and p'. pi are 
the LUKE grid points. The discrete form of fl^Tl) is then 

^_e'Ro^ n, jpi) , ^^^^ Ro [m] ^ jpi) [lO^'m-^] ^ 
^~4'Po^ ~ A[W] 4- r.(p.)[keV] ^^J' ^2^^ 

where /j is the current driven in the poloidal cross-section of Tj. Note that ( reflects 
the sign of the driven current. 

In Figures [6] -[9] we show the current drive efficiency ( for all the plasma and launch 
scenarios listed in Table [21 i.e., for the different frequencies, vertical launch positions 
and toroidal injection directions. The classification of the current drive mechanism is 
performed automatically by calculating the average (absorbed power weighted) N\\ of the 
rays and subsequently comparing the LUKE-calculated current direction to Ohkawa and 
Fisch-Boozer current directions. In certain cases, this leads to ambiguous results, either 
because the rays have mixed signs of N\\ or they are absorbed at different harmonics. 
The results were obtained by AMR and LUKE coupled simulations with 1 MW incident 
power. We immediately notice the importance of the launch parameters as they strongly 
influence the location of the wave power deposition (which obviously coincides with the 
driven current location) and the current drive efficiency in a fixed plasma equilibrium. 
Clearly, by changing these parameters, we can select a specific scenario — on/off axis 
deposition at almost any p with high/low \(\. There is full flexibility in the direction 
of the driven current because of the (a)symmetry (1191) . EBWs are most flexible and 
efficient in driving current in NSTX plasmas, mainly because in this case the magnetic 
field is monotonic without any magnetic well in the edge region. \(\ = 0.4 can be reached 
at almost any radius in NSTX. Our current drive efficiencies are similar to experimental 
values from COMPASS-D [6] or Wendelstein 7-AS [8] as well as to numerical results 
obtained for MAST-U |50] or NSTX [51]. 

While similar efficiencies can be reached using EC X- and 0-modes in the central 
region, the X- and 0-mode current drive efficiency typically decreases with radius, 
particularly because of trapping effects (see, e.g., |19]), which is not the case with 
EBWs. The L-mode plasma parameters cause higher absolute current drive efficiency, 
i.e., higher ri/(. There exist several significantly higher efficiency second harmonic cases 
with p = 0.1 (i.e., almost on the magnetic axis) in both L- and H-modes. Even though 
the results reported here indicate a limited on-axis accessibility and flexibility, this may 
not necessarily hold for slightly different frequencies or launch positions. 

Typically we flnd that in the central plasma regions we drive a Fisch-Boozer 
current |52] while Ohkawa current [53] in edge regions. Higher harmonic absorption. 
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i.e., absorption on the n harmonic with ncOce > favours the Ohkawa mechanism. 
Typically we can distinguish three EBW efficient current drive regions: 

(i) Fisch-Boozer current drive with lower harmonic absorption predominantly near the 
centre. 

(ii) Ohkawa current drive with lower harmonic absorption predominantly near the edge. 

(iii) Ohkawa current drive with higher harmonic absorption predominantly between the 
plasma centre and edge (i.e., between the ffist and second regions). 

The location and size of these regions are very different in the investigated plasmas 
and these regions typically overlap. There also exist cases when the harmonics 
are overlapping, i.e., the wave is absorbed on two different EC harmonics. Quite 
interestingly, in these cases the current is still driven in one direction even though the 
resonant electrons have their v^j-es having different signs. This occurs because the lower 
harmonic absorption favours the Fisch-Boozer mechanism (for which f Hres -j < 0), while 
the higher harmonic absorption favours the Ohkawa mechanism (for which fyres ■ j > 0). 

In going from an NSTX L-mode to an NSTX H-mode to a MAST-U and then to 
an NHTX plasma, the external magnetic field increases together with the appearance 
of magnetic wells near the edge (these wells being caused by strong edge currents) and 
we observe a decrease in the current drive efficiency, as well a decrease in the flexibility 
of the EBW absorption and central plasma accessibility. 
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Figure 6. Current drive efBciency C(synibols) and ri/( conversion factor (dashed line) 
versus p, NSTX L-mode first (full symbols) and second (open symbols) harmonics, 
all frequencies and vertical launch positions as listed in Table [21 both positive and 
negative -/V||o, 1 MW incident power. Neither the vertical launch position nor the iV||o 
sign can be graphically distinguished in the figure. 
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Figure 7. Same as Figure [S]- but for the NSTX H-mode. 



0.4 



0.2 



C 0.0 



-0.2 



-0.4 






















® ® 






® 

(B 


ll 


® * 






. ® 

® o ® 






•tr 

□ 




■ * 




n/C-' 













0.2 



0.4 



0.6 



0.8 



0.25 



0.20 



5 

0.15<| 



0.10 



.05 



frequency colours 


IS.O GHz 




32.0 GHz 


19.0 GHz 




33.0 GHz 


20.0 GHz 




34.0 GHz 


21.0 GHz 




35.0 GHz 


22.0 GHz 




36.0 GHz 


23.0 GHz 




37.0 GHz 


24.0 GHz 




38.0 GHz 


25.0 GHz 




39.0 GHz 


2G.0 GHz 




40.0 GHz 


27.0 GHz 




41.0 GHz 


2S.0 GHz 




42.0 GHz 


29.0 GHz 






30.0 GHz 






« Fisch- 


Boozer 


• Ohkawa 




* Udetermined 


■ Harmonic overlap 


L nco^g <u absorption 


H nu^g >b) absorption 



Figure 8. Same as Figure [6]- but for the MAST-U plasma. 



4-2. A^ii and quasilinear effects 

In order to show the effects of N\\, we calculate its average value (not to be confused 
with the initial N\\), weighted by the absorbed power: 

E A'||AP,ay(iV||,p) 

(Nn) = '-^^^^ . (24) 

i 

Figure [TO] shows the current drive efficiency versus [(A/^ij^l. It is found that the two 
quantities are clearly uncorrelated. As a consequence of the short wavelength of EBWs 
(/cj_pe ~ 1)5 the resonant v±_ is low, irrespective of the value of N\\. The dominant factors 
in EBW CD efficiency are the iV||-spectrum (mixing of signs), harmonic overlapping 
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Figure 9. Same as Figure [S]- but for the NHTX plasma. 



(because of large |iV|||) and Fisch-Boozer versus Ohkawa effects. Figures [TT] and [T2l show 
the current drive efficiency versus the absolute and the relative N\\ variance: again, we 
find no clear correlation. There is only a weak (logarithmic) decrease with the relative 
variance, starting at ~0.1. Most of the cases have a rather narrow N\\ spectrum, with 
absolute variance < 0.2 and relative variance < 0.1. NSTX H-mode, MAST-U and 
NHTX results show very similar behaviour. 



0.7 
0.6 

0.5l ^* 
0.4 
0. 
0. 
0. 

H'. 



.6 T T 1^* harmonic 

■5^* ' T i i 2"" harmonic 

_ T T A I 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

l(N||>l 



3.5 



Figure 10. Current drive efficiency versus the magnitude of the mean A^||, for all 
NSTX L-mode cases, Pq = 1 MW. 



We now compare the effect of two different vertical launch positions for the NSTX 
L-mode plasma at frequency 17 GHz. The ray trajectories and the evolution of N\\ are 
plotted in Figure [131 Those rays launched close to the midplane propagate straight to 
the magnetic axis, and the central ray's does not change appreciably until the ray 
gets close to the resonance (around i? = 1 m). | A^y | now starts to increase exponentially, 
and the beam splits in two parts that propagate in opposite vertical directions. Finally, 
the rays are absorbed, having been split in approximately two halves with opposite signs 
of iV|| at the absorption location. This behaviour demonstrates the typical behaviour of 
midplane rays at frequencies where cold EC resonance surface is lying in the inboard 
half of the plasma. For off-midplane launch at 17 GHz, one sees in Figure [13] that \N\\\ 
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Figure 11. Current drive efficiency versus absolute variance of A^||, for all NSTX 
L-mode cases, Pq = 1 MW. 
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Figure 12. Current drive efficiency versus relative variance of iV||, for all NSTX 
L-mode cases, Pq — I MW. 



steadily increases and the waves are absorbed at an EC resonance that has been Doppler 
shifted. Since all the rays have the same sign in iVy signs, one achieves a high current 
drive efficiency. 




R [m] R [m] 

Figure 13. Ray trajectories and the evolution of A^|| for two NSTX L-mode cases at 
17 GHz with different vertical launch positions; m and 0.3 m. 

These two examples clearly demonstrate how the deposition location and the 
current drive efficiency can be controlled by the choice of the vertical launch position. 
In Figures [H] - [17] we see the resulting power deposition and driven current densities, 
plotted for launched powers from 0.25 MW to 4 MW. The power deposition proffie (and 
consequently the driven current profile) is rather narrow in the midplane launch case 
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(Figures [m and dS]), since there is a sharp resonance close to the cold EC resonance 
surface. The driven current profile is oscillating around zero, resulting in nearly zero 
net driven current. These oscillations are rather artificial, partly because of the beam 
discretization by individual rays, and partly because, in reality, such oscillations would 
most probably be smoothed out by radial transport. There is almost no dependence on 
the launched power because of the sharp resonance. 
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Figure 14. Power deposition radial profile, NSTX L-mode, 17 GHz, 0.0 m vertical 
launch position. 
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Figure 15. Driven current density radial profile, NSTX L-mode, 17 GHz, 0.0 m 
vertical launch position. 



However, in the 0.3 m vertical launch position case (Figures [T6] and [TTl) . both the 
power deposition profile and the current drive profile are much broader. This is due to 
the strongly Doppler-shifted absorption with a relatively large N\\ spectral width, as well 
as Doppler broadening effects [3l] . The current is driven in one direction as the sign of 
A^ll is identical for all the rays. For this case there is power deposition on overlapping 
EC harmonics. Moreover, this is one of the interesting cases mentioned in the previous 
section, in which the Fisch-Boozer current from the deposition on the lower harmonic 
and the Ohkawa current from the deposition on the higher harmonic are in the same 
direction. 

In this configuration, the power is deposited on suprathermal electrons [31] and 
EBW H&CD is therefore strongly affected by quasilinear effects. Quasilinear flattening 
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of the distribution function with increasing power levels yields a relative reduction of 
the absorbed power, resulting in a more inward deposition. 
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Figure 16. Power deposition radial profile, NSTX L-mode, 17 GHz, 0.3 m vertical 
launch position. 
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Figure 17. Driven current density radial profile, NSTX L-mode, 17 GHz, 0.3 m 
vertical launch position. 

An input power scan for several NSTX L-mode cases is presented in Figures [18] and 
[T9l We find that there is no general tendency of the current drive efficiency to increase 
or decrease with the input power. Cases exist with increasing, decreasing or invariant 
( dependence on input power. However, in most cases, increasing power leads to either 
lower or similar current drive efficiency. In Figures [18] and [19] we also show the current 



profile maximum radial location pj and its width aj, defined as 

Pj max = arg max \j (p) | , (25) 
pe[o,i] 

Pj-i/2 = min {p : \j (p)| = \j {pjm^y,) /2| A ^ p < p^max} , (26) 

Pj+i/2 = max {p : |j (p)| = |j (pj ^ax) /2| A Pj^ax < p ^ 1} , (27) 

o^i = Pi+i/2 - Pj-l/2- (28) 



In other words, cXj corresponds to the full width at half maximum if the current profile is 
considered single-peaked. Increasing power causes the wave absorption to occur further 
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along the direction of propagation, which can either be towards the axis if the absorption 
occurs on the outboard side or away from the axis in the opposite case. This is caused by 
the quasihnear flattening of the distribution function and consequently lower absorption 
rate. 
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Figure 18. Input power scan of current drive efficiency and radial location of the 
current peak for 17 GHz NSTX L-mode cases with positive initial -/V||. Line-plots 
with symbols (colour online) represent ^ while symbols with vertical error bars (colour 
online) represent the radial current location pj ,„ax and its width aj (upper and lower 
limits represent Pj±i). Each group of the cross symbols belongs to one antenna 
position, and a horizontal shift is employed to separate the lines visually. 




0.1^ . — ■ ^0.8 




0.0 0.1 0.2 0.3 0.4 0.5 
antenna Z position [m] 



Figure 19. Same as Figure UHl but for negative initial A^j|. 



4.3. The effect of Z^s 

So far we assumed Z^s = 2, which is a realistic experimental value. However, Z^s can 
vary and it is very important to show the effect on EBW performance. We employ here 
the cases from the previous section and rerun the simulations with Zes ranging from 1 
to 3. Figure |20] shows the effect of Z^s on the current drive efficiency and the position 
of the current peak. Z^s affects the electron-ion collision frequency and particularly 
pitch-angle scattering. A larger value of Zes results in faster isotropization of current- 
carrying fast electrons. Thus, the current drive efficiency is inversely proportional to 
Zes [52] ■ The general trend of C versus Zee is shown in Figure EH Compared to Z^s = 2 
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results, the EBW current drive efficiency increases on average by 29 % for Z^s = 1 while 
a decrease of 18 % is observed for Z^s = 3. There is also a minor effect of Z^s on the 
EBW deposition location, as can be seen in Figure |20l This is again caused by the 
collision frequency change, which affects the plasma quasi-linear response to the wave 
power. 
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Figure 20. ZcB scan of current drive efficiency and radial location of the current 
peak for 17 GHz NSTX L-mode cases with negative initial iV||, Pq = 1 MW. Line- 
plots with symbols (colour online) represent C while symbols with vertical error bars 
(colour online) represent the radial current location pj max and its width aj (similarly 
to Figure [HI). 




Figure 21. Relative (with respect to Zcs — 2 values, denoted Co) changes of the 
current drive efficiency versus Zcs for all 17 GHz NSTX L-mode cases with both 
positive and negative initial A^|| , Pq = I MW. 



4.4- Robustness 

An important factor of any H&CD system is its robustness — the sensitivity to changes 
in plasma conditions or in the system's parameters. In the previous section we have 
already investigated what happens if the injected power is changed. Moreover, the effect 
of changing the vertical launch position can be seen in Figures [18] and [IHl this effect is 
rather strong and therefore the vertical launch position must be carefully chosen and 
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controlled. In this section we focus on the EBW H&CD performance sensitivity to 
plasma parameters. 

In Figures [22] and [23] we show the sensitivity of EBWs to plasma electron density 
and temperature variations in ±50 % range. All vertical launch positions and both 
initial N\\ signs of 17 GHz NSTX L-mode cases are used to calculate the medians 
of absolute location difference Apj^ax and relative current drive efficiency and profile 
widths |AC| / iCol and |Aaj| / |crjo|) where the subscripts denote results with the original 
plasma profiles. We first see a monotonic dependence of all the plotted quantities (except 
for two cases in |Acrj|), indicating a non-chaotic behaviour of EBW performance with 
changing plasma profiles. Quantitatively, the radial current location changes fractionally 
compared to the typical aj ~ 0.1. However, very precise localization might be important 
for certain applications, in which case a feedback system is highly advisable. The median 
difference in current drive efficiency is below 5 % for less than 25 % changes in the plasma 
profiles, which is very favourable. The current profile width is slightly more sensitive, 
a consequence of Doppler broadening. Not shown here are the variances. However, 
highest sensitivity is generally observed at lower frequencies, close to a midplane launch 
where the rays tend to oscillate, leading to current drive efficiencies that are typically 
low. In most cases the results are close to the median values. 
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Figure 22. Medians of absolute current location difference Apjniax and relative 
current drive efficiency and current profile width differences versus varying plasma 
electron density. 17 GHz NSTX L-mode 1 MW cases (similarly to Figure fT8)) are used 
to calculate the medians. 



Another parameter that can vary in tokamaks is the plasma current and the toroidal 
magnetic field. Unlike density and temperature profiles, which are only crudely pre- 
programmed and evolve during the discharge, it is typical that the plasma current 
and toroidal magnetic field do not change during the discharge (except, of course, in the 
start-up and shut-down phases) and that their properties are pre-programmed with high 
confidence. This makes the demands on the sensitivity on these quantities less stringent 
as compared to the temperature and the density. In Figures [2l] and [25] we show the 
sensitivity of 17 GHz L-mode cases to poloidal and toroidal magnetic field changes. The 
fields are simply changed by multiplying the respective components so that the resulting 
equilibrium is no longer a solution of the Grad-Shafranov equation. Significantly larger 
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Figure 23. Same as Figure [221 but for varying electron temperature 

effects of the magnetic field changes on EBW results can immediately be noticed. The 
sensitivity is particularly high for the toroidal field simply because the toroidal field is 
much larger than the poloidal field in most of the plasma cross-section. Also notice that 
changing the total magnetic field by 10 % is similar to changing the heating frequency 
by 1.4 GHz, which is the change in the central Ucc- For large magnetic field changes we 
can even change the EC absorption harmonic number — e.g., decreasing Stor by 25 % 
shifts 17 GHz into the second harmonic range. 
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Figure 24. Same as Figure [24j but for varying poloidal magnetic field. 



5. Conclusions 



By means of coupled ray-tracing and Fokker-Planck simulations, we have thoroughly 
investigated electron Bernstein wave H&CD prospects for spherical tokamaks. For the 
first time, a simple analytic formula for the 0-X conversion efficiency of a Gaussian 
beam is derived from ID plane wave theory. This formula supports our choice of the 
Rayleigh range as the antenna beam principal parameter that is fixed for all simulated 
cases. On an extensive set of EBW launch scenarios with varying frequency, vertical 
antenna position and toroidal injection angle, we show that EBWs can be absorbed at 
almost arbitrary radius and that EBWs can drive current with efficiencies comparable to 
electron cyclotron O- or X-modes. Moreover, the efficiency does not change with radius. 
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Figure 25. Same as Figure [24l but for varying toroidal magnetic field. Results for 
-Btor -50 % could not be calculated. 



while typically the efficiency of X- and 0-modes decreases with radius. Best results in 
terms of efficiency and flexibility are achieved in NSTX plasmas, where the electron 
cyclotron frequency radial profiles are monotonic. In general, normalized current drive 
efficiencies \(\ on the order of 0.3 - 0.4 are feasible for all target plasmas, absolute 
efficiencies then depend on the plasma parameters as /rfZ-Pq — O.SlCTg/i^o'^e, where 
the units are keV for T^, m for Rq and lO^^m"^ for rie- These results, however, do not 
reflect the mode conversion efficiency, which, as we have shown, is limited by the beam 
divergence and can be further degraded be other effects. 

For EBWs, the initial value of |A^||| is fixed by the mode-conversion process and 
only the sign of N\\ can be chosen at will. The evolution of N\\ is determined by 
the wave frequency, the vertical launch position and by the plasma parameters. We 
have shown how different vertical launch positions strongly influence the A'"]! spectrum 
and consequently the current drive efficiency. However, there seems to be no general 
correlation between the current drive efficiency and the iVy spectrum and its width. 

Input power scans have been performed to investigate the quasilinear effects. 
Increasing power generally leads to either lower or similar current drive efficiency. 
Higher power also causes the wave absorption to occur further along the direction 
of propagation, which can either be towards the axis if the absorption occurs on the 
outboard side or away from the axis in the opposite case. An important factor is the 
effective ion charge, which affects the electron-ion collisionality, and, consequently, the 
current drive efficiency significantly depends on Zgff . A minor effect of Z^s on the driven 
current location can be observed, which is caused by changing the plasma quasilinear 
response. 

The sensitivity of EBW H&CD to changes in plasma parameters has been 
investigated. It has been shown that the EBW performance is rather robust. Neither the 
current drive efficiency nor the radial location changes significantly when the electron 
temperature or density changes moderately. Larger sensitivity is observed for magnetic 
field changes, especially the (dominant) toroidal field. 

In conclusion, the EBW is a promising candidate for a powerful and flexible 
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auxiliary H&CD system for spherical tokamaks, in many aspects comparable to EC 
systems for standard aspect ratio tokamaks. However, several problems, which are not 
including in this survey, particularly the coupling, need to be addressed and better 
understood. High-power experiments on major spherical tokamaks or similar machines 
would be greatly beneficial to tackle these topics and validate our results. 
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